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Here we report a detailed study on spectroscopy, structure and dynamics of water molecules at 
air/water interface, investigated with Sum Frequency Generation Vibrational Spectroscopy (SFG- 
VS). Quantitative polarization and experimental configuration analysis of the SFG data in differ- 
ent polarizations with four sets of experimental configurations can shed new lights on our present 
understanding of the air/water interface. Firstly, we concluded that the motion of the interfacial 
water molecules can only be in a limited angular range, instead rapidly varying over a broad angular 
range in the vibrational relaxation time suggested previously. Secondly, because different vibrational 
modes of different molecular species at the interface has different symmetry properties, polarization 
and symmetry analysis of the SFG-VS spectral features can help assignment of the SFG-VS spectra 
peaks to different interfacial species. These analysis concluded that the narrow 3693cm~^ and broad 
3550cm~^ peaks belong to Coav symmetry, while the broad 3250cm~^ and 3450cm~^ peaks belong 
to the symmetric stretching modes with C2v symmetry. Thus, the 3693cm~^ peak is assigned to 
the free OH, the 3550cm~^ peak is assigned to the single hydrogen bonded OH stretching mode, 
and the 3250cm~^ and 3450cm~^ peaks are assigned to interfacial water molecules as two hydrogen 
donors for hydrogen bonding (with C2v symmetry) , respectively. Thirdly, analysis of the SFG-VS 
spectra concluded that the singly hydrogen bonded water molecules at the air/water interface have 
their dipole vector direct almost parallel to the interface, and is with a very narrow orientational 
distribution. The doubly hydrogen bond donor water molecules have their dipole vector point away 
from the liquid phase. 



I. INTRODUCTION 



Interfaces of water are the most important subjects not 
only because water is widely involved in physical, chem- 
ical, environmental as well as biological processes, but 
also because water is so far the most mysteries molecule 
in the universeiS*^ Among them, air/water interface 
has been intensively investigated theoretically or experi- 
mentally over the last decades. Spectroscopy, molecular 
structure and dynamics at air/ water interface is stud- 
ied with theoretical analysis such as ah initio calculation 
or molecular dynamics simulation^^iSiL&^iiSiiiiiSii^ or ex- 
perimental techniques such as X-ray reflection^^*i? Stim- 
ulated Raman Scattering (SRS)fi& Near-edge X-Ray Ad- 
sorption Fine Structure (NEXAFS)fii Second Harmonic 
Generation (SHG)fi24iS as well as Sum Frequency Gener- 
ation, etc.^°i^-^i^^i^'^i^^i^^i^^i^''i^^ Among these experimen- 
tal techniques, Second Harmonic Generation and Sum 
Frequency Generation are the most important methods 
for molecular interface studies because of their surface 
sensitivity and specificity. With these 
investigations, the properties of the water molecules at 
the interface, such as the surface density, surface struc- 
ture, surface potential as well as surface dynamics, have 
been intensively discussed. 
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However, conclusions on the surface molecule species 
at air/ water interface are still under discussion>i2i22i22i2& 
With SFG-VS experimental studies, the following inter- 
facial water species have been reported in literatures, 
namely, water molecules straddle at the interface with 
one OH bond hydrogen bonded to neighboring molecules 
in liquid phase (singly bonded OH) and another OH 
bond free from hydrogen bonding (free OH) in gas 
p]-^g^gg^22j23j24j37 -^^ater molecules with both OH bonds 
symmetrically hydrogen bonded in a tetrahedral network 
(ice- like and liquid- like structures) i^SiSiSiSl and water 
molecules in gas phase with both OH bonds not hy- 
drogen bonded pointing into the liquid phase^^'^^ With 
NEXAFS measurement and ab initio molecular dynam- 
ics simulation, water molecules with both OH bonds 
not hydrogen bonded pointing out of the interface was 
also proposediiiii The latter case is particularly contro- 
versial because NEXAFS is not strictly a surface spe- 
cific technique.^^ With polarization SFG-VS measure- 
ment, Wei et al. discussed the absence of SFG spectra 
in some polarization combinations and proposed an ex- 
planation through fast orientational motion in a broad 
range of about 102° in a time scale comparable or less 
than 0.5 ps^ However, puzzle still remains because 
some of the experimental studies suggests ordered and 
slow dynamics for interfaces of hydrogen bonding liq- 
uids, while some experimental investigations suggested 
a more dynamic and less ordered picture for the liq- 
uid interfaces, air/water interface includedi^ In addi- 
tion, whether the surface orientation relaxation is fast 
or slow than the bulk water molecules is also an issue 
under discussion in the recent literatures>iSi2Si^ Besides 
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SHG and SFG-VS experimental studies,^^ Structure 
and dynamics of water molecules at the air/ water inter- 
face have also been intensively discussed with theoretical 
simulationsi^iSiLSiSiiSiiiiiSii^ Even though with so much 
efforts and progresses both by experimentalists and the- 
oreticians, our detailed understanding of air/ water inter- 
face is still limited. Just as indicated by B. C. Garrett 
recently)^ '...(direct) experiments are difficult to perform 
because the liquid interface is disordered, dynamic, and 
small (typically only a few molecules wide) relative to the 
bulk'. 

Actually, direct measurement of the liquid interface 
is not as difficult as suggested as above. It has been 
known that along with SHG, SFG-VS can provide direct 
measurement on liquid interface no other technique can 
match j.^M.9r As pointed out by Miranda and Shen, 'SPG is 
currently the only technique that can yield a vibrational 
spectrum for a neat liquid interface'j^i In fact, with the 
advances of ultrafast laser and detection technology in 
the past decade and especially recent few years 
particularly with commercial systems designed for SFG- 
VS measurement,^^. SFG-VS, as well as SHG, experi- 
ments have come from easier to routine^i The real dif- 
ficulty lies on the fact that quantitative analysis and 
interpretation of the SFG-VS, as well as SHG, data 
had been not as well developed and widely performed 
until recently '5,35,48,49,50,51,52^53,54,55,56,57,58 Therefore, 

conclusions in many previous reports on the investiga- 
tions of air/water interface, as well as other liquid inter- 
faces, with SFG-VS are subjected to different interpreta- 
tions. 

As we have demonstrated in a series of re- 
cent publications, systematically quantitative treat- 
ment to SFG-VS data is not only possible, but 
also very effective for obtaining detailed spectroscopic, 
structural and thermodynamic properties of liquid 
interfaces.'^^i^Oi'^^i^'i^^i^^i'^^i^^i^^i^^ In these works, we not 
only developed methodology for quantitative polarization 
and experimental configuration analysis in SFG-VS and 
SHG, we also tested accuracy and sensitivity of some of 
the methodology. We have applied them to elucidated 
the anti-parallel double layered structure and thermo- 
dynamics of some organic liquid aqueous solution in- 
terfaces. In addition, we also demonstrated that a set 
of polarization selection rules (or guidelines) in SFG- 
VS can be developed for vibrational spectrum assign- 
ment through symmetry analysis of the SFG-VS spec- 
tral featuresi^Si^ This latter approach is extremely useful 
for discerning complex SFG-VS spectrum with unidenti- 
fied or controversial assignments. Recently, based on po- 
larization analysis, Ostroverkhov et al. demonstrated a 
phase-sensitive interference analysis of SFG polarization 
spectra of water/quartz interfaceiS With these develop- 
ment, in this report we intend to apply these analysis 
methodologies to the study of air/water interface. 

In this work, we examined SFG-VS spectra at 
air/ water interface measured in different polarizations 
under four experimental configurations with polarization 



analysis method and experimental configuration analy- 
sis. With these analysis, detailed new information are 
obtained for understanding of the spectroscopy, struc- 
ture and dynamics of the air/water interface. In the 
following sections, after a brief introduction of the theo- 
retical background and experimental conditions, we first 
discuss the motion of the interfacial water molecules at 
the air/water interface, which was previously suggested 
experiencing rapidly motion over a broad angular range 
in the vibrational relaxation time; then we use polar- 
ization and symmetry analysis of the SFG-VS spectral 
features for assignment of the SFG-VS spectra peaks; in 
the end, we shall discuss the structure and orientation of 
the water molecules at the air/water interface. 



II. POLARIZATION AND EXPERIMENTAL 
CONFIGURATION ANALYSIS IN SFG-VS 

Quantitative polarization analysis and experimental 
configuration analysis can provide rich and detailed infor- 
mation of spectroscopy, structure and dynamics of molec- 
ular interfaces.SSiS*^2iffi*aiS Generally, the SFG inten- 
sity in the reflective direction isj^^i^ 



I{co) = 



■K'^'ui sec (3 



c^ni{uj)ni{uJi)ni{LJ2) 



(2) 



IiLUl)I{LU2) (1) 



in which ui, coi and lu2 are the frequencies of the SFG 
signal, visible and IR laser beam, respectively. ni{uji) is 
the refractive index of bulk medium i at frequency t^i, 
and n'{uji) is the effective refractive index of the inter- 
face layer at cui. Pi is the incident or reflection angle 
from interface normal of the ith light beams; /(wi) is the 

(2) 

intensity of the SFG signal or the input laser beam. Xe/f 
is the effective second order susceptibility for an inter- 
face. The notations and the experiment geometry have 
been described in detail previouslyi '^^i^" 

(2) 

Xeff fo'^ fo^'^ generally used independent polar- 
ization combinations can be deduced from the 7 nonzero 
macroscopic susceptibility tensors for an achiral rotation- 
ally isotropic interface (Cooi>)i^^^ 

^(2y*p = Lyy{uj)LyyiuJi)L^,X^2)sinP2Xyyz 
^(2yps ^ Lyy{u;)L;,^{u;i)Lyy{ui2)sinPiXyzy 

^(2)^pss ^ Lzz{uj)Lyy{uJi)Lyy{uJ2)sinl3Xzyy 

x'eff^^ = -Lx^{uj)Lxx{LLii)Lzz{uJ2)cos(3cosf3isinP2Xxxz 
~Lxx {'-^)Lzz {'-^l)Lxx {uj2)cos(3sin(3icosP2Xxzx 
+Lzz{uj)Lxx{(^i)Lxx{(^2)sinf3cosl3icosl32Xzxx 
-\-Lzz{uj)Lzz{uJi)Lzz{uJ2)sinl3sinPisinf32Xzzz 

(2) 

It is so defined that the xy plane in the laboratory co- 
ordinates system X{x,y,z) is the plane of interface; all 
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the light beams propagate in the xz plane; p denotes the 
polarization of the optical field in the xz plane, with z 
as the surface normal, while s the polarization perpen- 
dicular to the xz plane. The consecutive superscript, 
such as ssp, represents the following polarization combi- 
nations: SFG signal s polarized, visible beam s polarized, 
IR beam p polarized, and so forth. La (i = x, y, z) is the 
Fresnel coefficient determined by the refractive indexes of 
the two bulk phase and the interface layer, and the inci- 
dent and reflected anglesi '^^i^° Xijk tensors are related to 

(2) 

the microscopic hyperpolarizability tensor (ilyj./ of the 
molecules in the molecular coordinates system A'(a, 6, c) 
through the ensemble average over all possible molecular 
orientations 



(2) 



i'j'k' 



{Rii'Rjj'Rkk')Piflk' 



(3) 



where Rw {0, 0, ■0) is the matrix element of the Eu- 
ler rotational transformation matrix from the molecu- 
lar coordination A' (a, 6, c) to the laboratory coordination 

\{x,y,z); Plijii^i is the microscopic (molecular) hyperpo- 
larizability tensoriSiSiSSi Here Ng is the molecular num- 
ber density at the interface. (A) represents orientational 
average of property A{9, 0, ^p) over the orientational dis- 
tribution function f(6, (/>, (p). 

/ ^^ _ lo C C MO, ^, i')f{0. 0, ^) sin 9d9 d^d^ 
/; /(^' ^) sin OdOdcj^d^ 

For SFG-VS, Z?*^^' is IR frequency {uj2) dependent. 



^(2) _ g{2) 
l-'i'j'k' ~ l~'NR,i'j'k' 



(2) 

Thus, Xijk "-^n be expressed into, 



E 



UJ2- OJq^- 



(5) 
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ijk 
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Xq,ijk 
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Therefore, SFG-VS measures the vibrational spec- 
troscopy of molecular interfaces. For dielectric interfaces, 

(2) 

such as liquid interfaces, the non-resonant term f3j^j^ ^i^iy 
is generally negligible compare with the reso- 



,(2) 



or XNR,.jk 
nant terms. 

Recently, we have found that the following formulation 
is very effective in quantitative polarization and orienta- 
tion analysis of SFG and SHG data. It can be generally 
shown that in surface SFG and SHG for an interface with 
orientational order, the effective second order susceptibil- 
ity Xe// '^^n be simplified into the following form4S 



X^^^f ^Ns*d* ((cos e)-c* (cos^ B))^Ns*d* r{9) (7) 



r{0) is called the orientational field functional, which con- 
tains all molecular orientational information at a given 
SFG experimental configuration; while the dimensionless 
parameter c is called the general orientational parameter, 
which determines the orientational response r{0) to the 
molecular orientation angle 9; and d is the susceptibility 
strength factor, which is a constant in a certain experi- 
mental polarization configuration with a given molecular 
system. The d and c values are both functions of the re- 
lated Fresnel coefficients including the refractive index of 
the interface and the bulk phases, and the experimental 
geometry. 

The key for quantitative analysis is that both d and 
c can be explicitly derived from the expressions of the 

(2) 

Xgff in relationship to the macroscopic susceptibility and 
microscopic (molecular) hyperpolarizability tensors for 
a particular molecular vibrational modesj^S*^ as shown 
for the water molecules with C2v symmetry in the ap- 
pendix. With the parameters c and d, the polariza- 
tion dependence and the orientation dependence of the 
SFG/SHG signal for a certain interface at certain ex- 
perimental configuration can be analyzed and calculated 
with clear physical picture on molecular orientation and 
orientational distribution!^ Reciprocally, information on 
the molecular symmetry, molecular orientation and dy- 
namics can be obtained from the analysis on the SFG 
intensity relationships measured in different polarization 
combinations and experimental configurationsi^SiSii^SiS 
The orientational average in Eq|31 is only the static 
average on molecular orientations, without considering 
fast molecular motion effects. Recently Wei et al. dis- 
cussed the fast and slow limit of the time average over 

(2) 

orientational motion for xlff ' ^nd they also applied this 
treatment to analysis the polarization dependence of SFG 
measurement of the OH stretching vibrational spectra for 
the air/water interface. ■^^ In the fast motion limit, the 
orientational motion is faster than the vibrational relax- 
ation time scale l/F^ of the qth vibrational mode; while 
in the slow motion limit, the orientational motion is much 
slower than l/Tq. 

According to Wei et al.,'^^ the slow motion limit gives. 



X 



(2) 



i = ^EE 



Pq,l'j'k' 



q i'j'k 

while the fast motion gives 



f-f UJ2-UJq+ iV, 



■{Rii'Rjj'Rkk') (8) 



xgi = ^EE 



^ qA'j'k' 



— 0J2 

q i'j'k' 



■ {Rw Rjj' ) {Rkk' ) (9) 



in which R\x'{t) — A • X'{t) is the time-dependent di- 
rection Euler transformation matrix from A'(a, 6, c) to 
X{x,y,z) coordinates system. Because of the molecular 
orientational motion, the molecular coordinates A' (a, 6, c) 
is time-dependent. Eq|S|is equivalent to Eq^l which is 
obtained by insertion of Eq[Sl into Eq|31 
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III. EXPERIMENT 

The details of the laser system has been described in 
our previous reportsi^Si^SiS, Briefly, the lOHz and 23 pi- 
cosecond SFG spectrometer laser system (EKSPLA) is 
in a co-propagating configuration. The efficiency of the 
detection system has been improved for the weak SFG 
signal of air/water interface. A high-gain low-noise pho- 
tomultiplier (Hamamatsu, PMT-R585) and a two chan- 
nel Boxcar average system (Stanford Research Systems) 
are integrated into the EKSPLA system. The voltage 
of R585 was 1300V in the measurement for air/water 
interface, and 900V for the Z-cut quartz surface. The 
wavelength of the visible is fixed at 532nm and the full 
range of the IR tunability is lOOOcm^^ to 4300cm~^. 
The specified spectral resolution of this SFG spectrome- 
ter is < 6cm~^ in the whole IR range, and about 2cm^^ 
around 3000cm~^. Each scan was with a 5cm~^ incre- 
ment and was averaged over 300 laser pulses per point. 
Each spectrum has been repeated for at least several 
times. Moreover, for sps polarization, each spectrum has 
been repeated for more than a dozen times and averaged. 
The energy of visible beam is typically less than 300/iJ 
and that of IR beam less than 150/iJ around 3000cto~^ 
and 3700cTO^^, and less than 100/iJ in the region in be- 
tween. These are comparable to literature reported val- 
ues for measurement of air/water interface.^" All mea- 
surements were carried out at controlled room temper- 
ature (22.0 ± 0.5°C) and humidity (40%) . The sample 
used was ultrapure water from standard Millipore treat- 
ment (18.2 MSI • cm). The whole experimental setup on 
the optical table was covered in a plastic housing to re- 
duce the air flow. No detectable evaporation effect was 
observed for SFG spectrum during each scan. 

The normalization procedure of the SFG signal in dif- 
ferent experimental configurations need to be specifically 
discussed. The detail of the normalization procedure for 
a single experimental configuration was presented in Xing 
Wei's Ph.D. dissertation,?^ However, the difference of co- 
herent length and Fresnel factors with different incident 
angles in the quartz SFG signal measurements has to be 
corrected when comparing SFG signal in different exper- 
imental configurations. Therefore, the measured spec- 
trum is firstly normalized with the energy of the incident 
laser beams, and then normalized to the SFG signal of 
Z-cut quartz (also normalized by the energy of the in- 
cident lasers). Then it times with a converting factor 
between different experimental configurations. This fac- 
tor contains the influence of the coherent length of Z-cut 
quartz, the Fresnel coefficients,^'^ the Xijk value for Z- 
cut quartz, and the factor sec^P for each experimental 
configuration. Therefore, the end result is directly pro- 
portional to the SFG intensity in Eq^ If the spectrum in 
Fig. n]is divided by the factor and the factor of the 

PMT efficiency between 1300V and 900V, which is de- 
termined as 24.1 in our detection system, and then times 
the unit factor 1 x 10~'^^V'^m~^ which we left out for 
simplicity of graph presentation, it will give the value for 



|Xe//P- For example, the peak at about 3700cto ^ in the 
ssp spectra of Config.2 in Figgis about 0.23 unit. After 
above conversion it gives jXg^jp — 4.7 x 10^'^'^y^m^^, 
matching satisfactorily with the reported value for less 
than 10% differencei^ 

Even though the normalized intensities are generally 
consistent with each other, there can be possibly other 
sources of errors when intensities in different experimen- 
tal configurations need to be compared. Because the vis- 
ible and IR beams have different coherent lengthes in the 
Z-cut crystal, and because these coherent lengthes vary 
with different experimental incident angles, one of the 
most likely error might come from the different focus- 
ing parameters with different beam overlapping quality 
of the visible and IR beams in the Z-cut quartz crys- 
tal with different experimental configurations. There- 
fore, quantitative comparison of the SFG spectral inten- 
sities in different polarizations with the same experimen- 
tal configuration can be more accurate than comparison 
intensities between different experimental configurations. 
Even though the latter is a good solution to reduce such 
relative error associated with different experimental con- 
figurations need to be developed. 

IV. RESULTS AND DISCUSSION 

A. Polarization SFG Spectra of the air/water 
interface 

Firstly we would like to present the polarization SFG 
spectra of the air/water interface measured in four dif- 
ferent experimental configurations. 

We have demonstrated recently that the change of the 
SFG spectra in different polarizations by varying the ex- 
perimental configurations can be used for quantitative 
polarization analysis and orientational analysisi^SiS Here 
we present in Fig^the SFG spectra in the ssp, ppp and 
sps polarizations on the air/water interface at four ex- 
perimental configurations with different incident angles 
for the visible and IR laser beams. They are, Config.l: 
Visible=39°, IR=55°; Config.2: Visible=45°, IR=55°; 
Config.3: Visible=48°, IR=57°; Config.4: Visible=63°, 
IR=55°. 

There are four apparent peaks can be identified in 
the SFG spectra in Fig^ They are around 3700cm~^, 
3550cr7i~^, 3450cm~^ and 3250cto~^, respectively. The 
3700cTO~^, 3450cm~^ and 3250cto~^ peaks has been 
extensively discussed in the SFG literature iSiiSSiSiSiSi 
However, the 3550cm~^ peak has been observed, but not 
yet clearly identified or assignedi^^ The results of global 
fit of these spectra with four Lorentzian peaks in Eq|^ 
are listed in Table From the fitting results we can see 
that the peak bandwidths of the 3550cto"^, 3450cto~^ 
and 3250cm~^ peaks are 77 ± llcm~^, 103 ± 7cm~^ and 
89 ± 9cm~^, respectively. Such broad bandwidths indi- 
cate that they all belong to different hydrogen bonded O- 
H stretching vibrational modes. However, the bandwidth 
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3000 3200 3400 3600 3000 3200 3400 3600 3000 3200 3400 3600 3800 

Wavenumber (cm ) 



FIG. 1: SFG spectra of air/water interface in different polarization combination and experimental configurations. All spectra 
are normalized to the same scale. The solid lines are globally fitted curves with Lorentzian line shape function in Eq. |S| Note 
the different error bars for graphs in different scales. 



of the 3693cm~^ peak width is only 17cm~^, consistent 
with the symmetric stretching (ss) vibrational mode of 
the free 0-H bondi^ The signs in Table contain the 
information of the relative phase and interference effects 
of the different vibrational modes. Here the phase of the 
3693cr7i~^ peak is held positive in each fit. Altering the 
relative phases of the peaks on the same spectrum can 
not give a reasonable fit. Because we used global fitting 
with all the spectra, these relative phases can be deter- 
mined accurately. They can be used to determine the 
symmetry properties of each vibrational mode in Section 
IV. C. 

According to Eq[21 the ssp spectra in different experi- 
mental configurations should have the same features from 



the Xyyz term. As shown in Fig[21 all ssp curves over- 
lap quit well when normalized to the 3693cm~^ peak. 
Calculation of the Fresnel factors with different incident 
angles can quantitatively explain the relative intensities 
in all four configurations.^^ Because the SFG spectral in- 
tensity from the air/water interface in the OH region is 
usually several times smaller than that of the C-H region 
from other air/liquid interfaces, the air/water interface 
SFG spectra are usually very hard to measure experimen- 
tally. Therefore, the well overlapping of the ssp spectra 
in different experimental configurations is a proof for the 
quality of our SFG-VS data. Furthermore, the spectra we 
obtained agree very well with these in the literatures 

In principle, the sps spectra in different experimental 
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TABLE I: The fitting results of tlie SFG spectra at air/watrer 
interface. Tfie spectra are fitted witli Lorentzian line shape 
function as Eq|S| The peak position of the vibrational modes 
LVq, the peak width Tq and the oscillator strength factor 
Xeff,q,ijk of the vibrational modes are listed. The first col- 
umn is the fitted value for XNR,eff,ijk- The relative error in 
fitting of sps is larger because of the small signal strength for 
sps spectra. 
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configurations should also overlap with each other when 
normalized. However, consistent with the calculations 
of the corresponding Fresnel factors, the sps signal level 
for Config. 1, 2 and 3 are very close to the noise level, 
and features in the sps spectra can not be clearly iden- 
tified except for the spectra of Config. 4. Therefore, such 
normalization and comparison for sps spectra is not as 
meaningful as the ssp spectra. 

Different from the ssp and sps spectra, the features in 
the PPP spectra in Fig^ changed drastically with differ- 
ent experimental configurations. This is because that the 
PPP spectra is determined by combination of four differ- 
ent Xijk tensors. Detailed polarization analysis and ex- 
perimental configuration analysis of these changes in the 
PPP spectra can provide symmetry properties for each 
spectral features, as well as orientation and structure in- 
formation of the interfacial molecular groups, as shall be 
shown lateri^2i^iiS We shall show that analysis of the 
PPP spectra in different experimental configurations is 
very informative. However, this advantage of ppp spec- 
tra analysis has not been well utilized in the previous 
literatures. 



B. Orientation and Motion of the Free OH Bond 

Now with the knowledge of the SFG vibrational spec- 
tra of the air/water interface, we can discuss the orien- 
tation and motion of the free 0-H bond at the air/water 
interfaced 



O sspl 

0.25-1 □ SSp2 
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V ssp4 




3200 3400 3600 
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FIG. 2: Overlap of the normalized ssp spectra of the air/ water 
interface in different experimental configurations. 



The sharp peak around 3700cm~^ was generally ac- 
cepted as the free OH bond protruding out of the liquid 
.^g^|-gj.^20j23j25j66 g^j^^j j^^^g ^ggjj treated with Coov sym- 
metry in polarization analysis i22i2^ Wei et al. studied the 
polarization dependence of the intensity of this peak in 
the ssp, PPP, and sps polarizations measured with exper- 
imental configuration of Visible^ 45°, IR=57°iS^ Their 
SFG-VS data are quantitatively very close to our data 
with Config. 2 as expected. Therefore, the ssp intensity 
of the 3693cm~^ peak is about 10 times of that of ppp, 
and the sps intensity is essentially close zero. Wei et 
al. realized that using the step orientational distribution 
function in Ea llOl as well as other distribution functions, 
such as Gaussian, centered at the surface normal, can not 
explain such ssp, ppp and sps intensity relationships with 
the slow motion average in EqlHl On the other hand, the 
fast motion average centered at the interface normal, as 
shown in Eq^with 6m = 51°, can fairly well explain the 
observed intensity relationships. Thus, it was concluded 
that the orientation of the free OH bond of the inter- 
facial water molecule varies over a very broad angular 
range {9m = 51°) within the vibrational relaxation time 
as short as Q.bpsi^ 
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As shown in Fig|31 Wei et a/.'s treatment predicts 
clearly zero intensity for the sps spectra at 3693cto~^ 
with the assumption of fast orientational motion centered 
at the surface normal. Using exactly the same parame- 
ters, our calculation of Config. 2 gives the same results 
as that by Wei et al. as it should have heenm^ It is 
clear that the simulation results in FigOlcan fairly well 
explain the data in Config. 1, 2 and 3, because all of 
them have relatively very small sps spectral intensity at 
3693cr7i~^. However, even though the fast orientational 
motion picture can explain the relative intensity between 
the ssp and ppp polarization in Config. 4, it is clear that 
it can not explain the apparently non-zero sps intensity 
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FIG. 3: SFG intensity of the free OH bond simulated with 
both slow motion limit (solid curves) and fast motion limit 
(dotted curves) following the procedure and parameters as 
Wei et alw~. 9m is the range of orientational motion of the 
free OH bond. All the curves presented include the factor of 
sec^ /3, and all intensities are normalized to the ssp intensity 
in Config.4 with 6m — 0°. The vertical lines indicate the 
distribution width suggested by Wei et al. 



at 3693cm~^ with Config.4. As long as the orientation 
distribution or orientational motion is assumed to be cen- 
tered to the interface normalj^^ orientational distribution 
functions other than the step function in Eg 1101 give the 
same conclusion. Since the slow motion limit is already 
not an optioufSi alternative description of the motion and 
orientation at the air/water interface has to be invoked. 

Because the air/water interface is rotationally isotropic 
around the interface normal, now we assume that the 
molecular orientation is centered around the tilt angle 
6o 0, instead of the interface normal {Oq = 0). If the 
Gaussian distribution function is assumed, we have 



1(0) 



(11) 



TABLE II: The general orientational parameter c and the 
strength factor d for the vibrational stretching mode of free 
OH bond in different experimental configurations. The d 
value bear the unit f3ccc of single OH bond. 

d-ssp c-ssp d-sps c-sps d-ppp c-ppp 
Config.1 0.274 0.515 0.112 1 -0.154 1.53 
Config.2 0.256 0.515 0.118 1 -0.120 2.05 
Config.3 0.248 0.515 0.118 1 -0.104 2.43 
Config.4 0.176 0.515 0.107 1 -0.035 6.55 



The details of the calculation of c and d can be found 
elsewhere i^2i^i*S It is clear from Tablemthat the c values 
for the ssp and sps polarizations are the same for all 
four experimental configurations; whereas the c values 
of the ppp polarization differ significantly for different 
experimental configurations. 
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FIG. 4: The simulated SFG intensity of vibrational stretch 
mode for free OH bond at different orientation angle 6 as- 
suming (7 = 0°. The factor sec? (3 in EqQis also included for 
comparison of SFG intensity in different experimental config- 
urations. All curves are normalized to the ssp intensity in 
Config.4 with 6^0 = 0°. The vertical lines indicate the orien- 
tation which quantitatively explains the observed SFG data. 



in which a is the standard deviation of the angular dis- 
tribution. We shall show in the foUowings that by using 
this distribution function, the 3693cto~^ peak in differ- 
ent polarization and experimental configurations in Fig. 
^can be quantitatively analyzed. 

Because the 3693cm~^ peak belongs to the ss mode 
of the free O-H bond at the air/water interface, it has 
been treated with Coov symmetry. Now we calculate 
the general orientational parameter c and the strength 
factor d for the ssp, sps and ppp polarizations in all four 
experimental configurations with the same parameters 
of the air/water interface as those used by Wei et al?- 



As we have demonstrated previously;^2i^2i^i*^ the [d * 
r{9)]'^ vs. 9 plot with cr = in different polarizations 
can provide direct first look of the physical picture for 
polarization analysis of SFG-VS data. Here we plot 
[d * r{6) * sec PY vs. 9 in FigQlin order to compare data 
in different experimental configurations. Thus, the rel- 
ative intensity for the 3693cm^^ peak in experimental 
Config.1, 2, 3 and 4 can be used to calculate the ori- 
entation angle of the free O-H bond. Using the known 
procedure3^2iS2*^ii^ and parameters, they give the fol- 
lowing four values, i.e. 28.7±1.2°, 32.6±0.5°, 34.6±0.7° 
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and 35.8 ± 1.0°, respectively. These values agree with 
each other quite well. However, the value from Config.l, 
whose ppp and sps intensities are both very weak, is not 
as reliable as the other three configurations. Averaging 
over these values gives 9 = 33° ± 1°. 

It is clear that cr = 0° is not physically possible for the 
liquid interface. However, the apparent success of the 
quantitative explanation of the observed SFG spectra of 
the free 0-H bond in different experimental configura- 
tions using a — 0° indicates that the actual a value can 
not be very broad. Simulation of the 3693cr7i~^ peak in 
different polarizations in each of the four experimental 
configurations using the Gaussian distribution function 
in Eq^2 concludes that a has to be smaller than 15° 
to satisfy the measured 3693cm^^ peak intensities in all 
four experimental configurations, cr — 15° is the largest 
distribution width allowed by the SFG experiment data 
for a Gaussian orientational distribution function. With 
a = 15°, we have ~ 30°. This indeed confirms that 
the orientation of the free 0-H bond is within a rela- 
tively narrow range (between 30° to 33°), with a rela- 
tively small distribution width (cr < 15°). Calculation 
with both Eq|Sl i.e. slow average limit, and EqO i.e. 
fast average limit, gives indistinguishable results with a 
as small as < 15° if 9o is around 30°. This is because that 
with a small distribution width, fast and slow motion av- 
erage should be the same according to Eq|51 and Eq|Sl 
Using a step distribution function around 6*0 ^ 0° also 
give very close orientation angle and distribution width. 

Thus, our conclusion of the free 0-H orientation and 
distribution at the air/ water interface is drastically dif- 
ferent from the conclusion given by Wei et at, which 
concluded that the free 0-H bond orientation varies in 
a broad range as big as 102° and as fast as 0.5 picosec- 
ond, which is the relaxation time for the 0-H stretch- 
ing vibration. It is clear that our conclusion is based 
on the successful explanation of the observed polarized 
SFG spectral intensities in different experimental config- 
urations, especially the relatively small but clearly non- 
vanishing SFG spectral intensity at 3693cto~^ in the sps 
polarization. Our conclusion explicitly supports ultrafast 
libratory motions with a relatively narrow angular range. 
As we have known, the dynamics libratory motion of the 
hydrogen bonding can be as fast as 0.1 picosecond.-^'®* 
Even with such ultrafast dynamics, the air/water in- 
terface is nevertheless well ordered. This is consistent 
with the generally well ordered picture of the air/liquid 
and air/liquid mixture interfaces. Recent quantitative 
analysis of data in SFG vibrational spectroscopy have 
suggested that vapor/liquid interface are generally well 
ordered, and sometimes even with anti-parallel double- 
layered structure3Si5SiSLS2iSMa. 

It has been generally accepted that liquid interface 
with strong hydrogen bonding between molecules should 
be well ordered^i. Our analysis here not only confirmed 
this conclusion, but also provided solid and direct exper- 
imental measurement of the orientation and motion at 
the air/water interface. 



C. Polarization Analysis and Determination of 
Spectral Symmetry Property 

Here we try to apply polarization analysis for identi- 
fying the symmetry property and for assignment of the 
SFG vibrational spectra of the air/water interface. 

The assignment of the SFG spectra of air/water inter- 
face in the range of 3000 to 3800cto~^ has been discussed 
intensively,'^ii'^^.^3,24,26,27,28,37,7Q,7i72 Richmond recently 

reviewed the current understanding of the bonding and 
energetics, as well as the SFG spectra assignment, 
of various aqueous interfaces, including the air/water 
interfacei2L2& The SFG spectral assignments heavily re- 
lied on band fitting of IR and Raman peak positions of 
bulk water or water cluster spectra^S as well as based 
on theoretical calculationsi2i2iZi24 The sharp peak at 
about 3700cm~^ has been unanimously assigned to the 
free 0-H stretching vibration mode. The broad peaks 
around 3250cm~^ and 3450cto~^ undoubtedly belong to 
the hydrogen bonded 0-H stretching modes, but their 
assignments are not as unanimous as the 3700cto~^ 
peak. The spectrum around 3250cto~^ was assigned to 
a continuum of 0-H symmetric stretches(ss), vi of wa- 
ter molecules in a symmetric environment (ss-s) , and was 
generally referred as " ice- like" region because of its sim- 
ilarity in energy to 0-H bonds in bulk ice. The broad 
band around 3450cm^^ was assigned to more weakly 
correlated hydrogen bonded stretching modes, and was 
called the "liquid- like" hydrogen-bonded region, where 
water molecules reside in a more asymmetrically bonded 
(as) water environment The broad peak around 
3550cTO~^ appeared clearly in the ppp SFG spectra has 
been identified once and it has not been clearly assigned 
so fariS^ Shultz et al. pointed out that these broad 
peak should also include the asymmetric stretching mode 
of water molecules in a symmetry environment and the 
bending overtonei24 Richmond et al. also suggested that 
the intensity at about 3450cm^^ include the contribution 
of donor 0-H bond.^® 

Recent progresses on SFG-VS have made it pos- 
sible to determine the symmetry properties of SFG- 
VS vibrational spectral features through comparison of 
SFG spectra in different polarizations and experimen- 
tal configurations i^S'Si'SiS The key idea of this develop- 
ment is from the commonsense of molecular spectroscopy 
that vibrational modes of molecular groups with differ- 
ent symmetry properties have different polarization de- 
pendence on the interacting optical fields iT^'7^ Applying 
these ideas to polarization analysis of SFG spectra has 
led to a set of polarization selection rules for different 
stretching vibrational modes of molecular groups with 
different molecular symmetry properties, such as stretch- 
ing vibrational modes of the CH3 (Ca^), CH2 (C21,) and 
CH (Coov) groupsi^SiSiiS Many of these selection rules 
are independent from molecular orientation and orienta- 
tional distribution. Therefore, they can be directly used 
to identify symmetry property of SFG stretching vibra- 
tional band. These progresses make it possible to analyze 
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FIG. 5: The simulated SFG intensity for symmetric stretching (ss) mode (left) and asymmetric stretching (as) mode (right) of 
water molecule with C2v symmetry. All curves presented include the factor of sec^ /3. The intensity of ss mode is normalized to 
the ssp intensity in Config.4 with ^ = 0° . The intensity for as mode is normalized to the sps intensity in Config.4 with 60 — 0° . 
The units between plots of the ss and as modes differ by 9.11 times according to the Pccc and jSaca values in the appendix. 



SFG vibrational spectra in situ, instead of rely only on 
the assignments from Raman and IR studies of the bulk 
phases, which can be called ex situ. Because SFG spectra 
usually has more features than those from IR and Raman 
measurement, some confusions and errors in the previ- 
ous spectral assignments have also been clarifiedi^SiSii^ 
Even though SFG is naturally a polarized spectroscopy 
and the interfacial molecular groups are usually ordered, 

this idea has not been systematically explored until very 
recently|49,50,5i,52 

Water molecule possesses C2V symmetry. If the two O- 
H bond of a water molecule are asymmetrically bonded, 
both O-H bond has to be treated separately with Coov 
symmetry. This classification of the water molecule sym- 
metry is generally true no matter it is hydrogen-bonded 
or not, in cluster, in bulk or at the interface. Therefore, 
the symmetry property of the SFG vibrational spectra 
features of the air/water interface can all be classified 
accordingly. Thus, there are three kind of stretching vi- 
brational modes for us to deal with, namely, the sym- 
metric (ss) and asymmetric (as) stretching modes for C21, 
symmetry, and the stretching mode for Coov symmetry. 
It is fairly easy to distinguish these three stretching vi- 
brational modes from the polarization selection rules for 
SFG spectra of CH2 and CH groups i^Si^ii^ Because the 
bond angle of the CH2 group is slightly larger than that 
of water molecule, the polarization dependence of the C2V 
water molecule are slightly different. The key difference 
is that even when the water molecule at the interface 
rotate freely around its symmetry axis, the sps spectral 
intensity of its ss mode does not vanish as that for CH2 ■ 
However, this fact does not make the polarization selec- 
tion rules different for the interfacial CH2 group and the 
C2V water molecule. 

Two of the major selection rules for the C2V group at a 
dielectric interface are: (a) ssp intensity is always many 



times of that of ppp for ss mode, and (b) ppp intensity 
for as mode is always several times of that of ssp. That 
is to say, if there is any peak which is stronger in the 
ssp than ppp spectra, it can not be from the as mode^^^ 
These two rules are independent from molecular orien- 
tation and orientational distributions at a rotationally 
isotropic interface. 

It is clear from these selection rules, the sharp peak 
around 3693cm~^ in Fig^ does not belong to the C2V 
symmetry, because its intensity in ppp polarization in 
Config. 1,2,3 is smaller than that in ssp polarization; 
while larger in Config.4. On the other hand, these all fits 
well with the simulations in Fig^ Therefore, 3693cto~^ 
peak is with Coov symmetry as the free O-H stretch- 
ing mode. Dissimilarly, both the broad peaks around 
3250cto~^ and 3450cm~^ in Fig^are very strong only in 
the ssp spectra in all experimental configurations. They 
fit well with the ss mode of the C2V symmetry, and can 
not belong to the as mode of the C2V symmetry, or the 
Coov symmetry. 

It is not so easy to determine the symmetry property 
of the broad 3550cm~^ peak, because it appears to be 
buried in the high frequency tail of the broad 3450cto~^ 
peak. It is not so straight forward to read its relative 
intensity in ssp and ppp polarizations from Fig^ How- 
ever, it appears significantly bigger in sps polarization in 
Config.4 than that in Config. 1,2, 3. Therefore, it appears 
to fit with the simulations in Fig^ In order to exclude 
the possibility that it may belong to a C2V mode (ss or 
as) , detailed simulation of the C2V modes in different po- 
larizations and experimental configurations is now called 
upon. 

As describe in the appendix, the parameter c and d of 
the C2V vibrational stretching modes are calculated for 
different polarizations and experimental configurations. 
Plots of [d * r{6) * sec/3]^ vs. tilted angle 6 of the water 
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molecule c axis from the interface normal using these c 
and d values are presented in Fig|Sl These plots again 
confirm that the 3693cm~^ can not belong to any C21, 
mode, especially with the one order of magnitude increase 
of this peak in ppp polarization. 

Clearly, the polarization dependence of the broad 
3550cr7i~^ peak does not fit to the ss mode of the C2V 
symmetry. Otherwise, according to Fig|5l its ppp inten- 
sity in Config.3 has to be about one order of magnitude 
weaker than that in the observed spectra. This peak can 
not be the as mode of the C2V symmetry either. Accord- 
ing to the c and d values for the as mode of the C21, 
symmetry, the phases in ssp and ppp polarizations have 
to be with opposite signs in all four experimental con- 
figurations. However, fitting of the ssp and ppp spectra 
indicates that in Config.4, the oscillator strengths had the 
same signs for ssp and ppp polarizations, even though in 
Congfig. 1,2,3, the oscillator strengths of this peak do 
possess opposite phases these two polarizations. This in- 
dicates that the 3550cto~^ peak is not C2V symmetry, 
and it appears to have Coav symmetry. 

Because in the bulk phase there is no observation of 
free 0-H bond, and because this broad peak 3550cto^^ 
appears to be hydrogen-bonded, there is only one possi- 
bility that it is the other 0-H bond of the interfacial wa- 
ter molecule which has a free 0-H bond extruding away 
from the liquid bulk phase. According to Fig^ the Coov 
0-H stretching mode in sps polarization is about twice 
as large in CogfLg.4 as that in the other configurations. 
This is fully consistent with the SFG spectra data in 
FigHand the fitting results in Table Furthermore, in 
Table U the phase of the broad 3550cto~^ peak is just 
opposite to that of the 3693cto^^ peak in both ssp and 
sps polarizations, indicating these two 0-H pointing to 
opposite directions. The phase of the ppp polarization 
of the broad 3550cm^^ peak changes signs with different 
experimental configuration. This is because the orienta- 
tional angle of the two 0-H bonds are some times on the 
same side of the minimum on the ppp curves in Fig01 
and sometime on the different side of the minimum, just 
as predicted with the experimental configuration anal- 
ysis. These detail features indicated the ability to un- 
derstand very subtle dependence of the SFG spectra on 
experimental configurations and the parameters used for 
the spectra calculations. Further study shall be reported 
elsewhere. 

Further support for the assignment of the broad 
3550cTO~^ peak came from the IR spectra measure- 
ment of the water dimer clusters, where the stretch- 
ing frequency for the donor 0-H bond is just at about 
3550cr7i~^iSSi2Si2L2& This assignment is a good support 
for our assignment of the peak at 3550cm~^ in ppp spec- 
tra to the single hydrogen-bonded water molecule at 
the interface. Furthermore, the two 0-H stretching vi- 
brations for the methanol dimer are at 3574cto~^ and 
3684cm^^i2^ The donor 0-H stretching mode is also in 
the same region of 3550cm~^. 

There is no observable spectra features in Fig^for the 



as mode of the C2V water molecules, neither hydrogen 
bonded nor non-hydrogen bonded. According to Fig|31 
for the as modes corresponding to the ss mode around 
3250cr7i^^ and 3450cto~^, their intensities have to be at 
least one order of magnitude weaker than that of their 
corresponding ss modes. It is understandable that we do 
not observe them. Above discussion also throw doubts 
on the existence of interfacial water molecules with two 
free 0-H bonds, as suggested somewhat less convincingly 
by some recent studies li^ii^iSiSSiii According to the po- 
larization selection rules and the calculation for the po- 
larization dependence of the C2V water molecules, no de- 
tectable spectral features satisfying the C2V symmetry in 
the 3600cTO~^ and 3800cm~^ has been observed in the 
SFG spectra. 

Here we clearly see that how polarization selection 
rules, quantitative polarization and experimental config- 
uration analysis can help determine the symmetry prop- 
erty of the observed spectra features. The importance 
of studying of spectral interference has been demon- 
strated in recent reportsi^Si^MS Analysis in this work also 
demonstrated that, in order to discern spectral details, 
it is useful and effective to analyze the spectral interfer- 
ence of different spectral features through global fitting 
of SFG spectra in different polarizations and experimen- 
tal configurations, and to compare fitting results with the 
prediction from the calculated c and d values. This also 
indicates the usefulness of the formulation of total SFG 
signal with functions of c and d parameters in Eq|7| 



D. Molecular Structure at Air /Water Interface 

With the analysis of the orientation and motion, vibra- 
tional spectral symmetry of the water molecules at the 
air/ water interface in above sections, we can have more 
understanding of the molecular structure of the air/water 
interface. 

In Section IV. B, we have determined that the free O- 
H oriented around 30° away from the interface normal 
with a orientational distribution narrower than a = 15°, 
and in Section IV. C, we have identified the spectral fea- 
ture around 3550cr7i^^ of the other hydrogen bonded 0-H 
bond of this water molecule. If the plane of this interfa- 
cial water molecule is close to perpendicular to the inter- 
face, the orientation of the hydrogen-bonded 0-H should 
point into the liquid phase with a orientation around 
135° away from the interface normal. This orientation is 
fully consistent with the calculation of the polarization 
and experimental configuration dependence of the broad 
3550cr7i~^ peak with a Coov symmetry with the observed 
SFG intensities, detail to be reported elsewhere. Such 
orientation makes the dipole of this water molecule points 
around 97° from the interface normal. This picture is 
fully consistent with conclusions in many previous experi- 
mental and theoretical studiesr''''^'^i'°''^i^''''^i^^i^^i^^i^^i^^ 
but certainly different from some^ 

From Section IV. C, the broad spectral features be- 
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tween 3100cto~^ to 3500cto~ are determined to be sym- 
metric stretching modes of the symmetry. Because 
the peaks are broad, and their energies is in the range of 
hydrogen-bonded 0-H stretching range, they can only 
come from the water molecules with two donor 0-H 
bonds, whose oxygen atom can accept either two, one 
or zero hydrogen atom from other water molecules as hy- 
drogen donors. Certainly, the water molecule with the 
oxygen atom forming two hydrogen- bonds is tetrahedral 
in shape and is "ice- like". This is consistent with the 
previous assignment of the broad 3250cm"^ peak. 

The water molecule with no hydrogen bond for the 
oxygen atom is obviously with C2V symmetry. How- 
ever, the water molecule with only one hydrogen bond 
for the oxygen atom may or may not preserve the 
symmetry. However, if this hydrogen bond perturbation 
to the water structure is limited, this water molecule 
can still be treated as with C2V symmetry. The last 
two kinds of water molecules are certainly not "ice- 
like" , but " liquid-like" . The two " liquid-like" species 
may have slightly different 0-H vibrational frequen- 
cies. However, only two apparently broad peaks in the 
BlOOcTO"^ to 3500cr7i~^ region have been identified in the 
literaturesi^^i^^i^'^i^4i^6,27,28,37,70,7i,72 ^ore studies on the 

possible hydrogen-bonded species are certainly warranted 
in the future. 

Here we confirm the conclusion by Brown et al. that 
these C2V water species all have their dipole vector point 
out of the bulk liquid phase, i.e. with both hydrogen 
atoms point into the bulk liquid phase. It is clear in Ta- 
ble U the signs of the ssp polarization oscillator strength 
factors of the C2V water species are all in opposite phase 
to that of the free 0-H peak at 3693cm~^ in all experi- 
mental configurations. The signs and values of the c and 
d parameters of the C2V and Coov in Table Hill and Table 
IhI respectively, indicate that the c axis of the C2V species 
has to be in opposite direction to the c axis of the free 
0-H bond at the interface. Therefore, as defined as in 
the appendix, the C2V species have to have their dipole 
vector point out of the bulk liquid phase. The calculation 
of the phase of the ppp as well as sps spectral features 
are all consistent with this picture. 

However, because the SFG spectral intensities of the 
ppp and sps polarizations are generally in the noise level 
in the 3100cm^^ to 3500cr7i~^ region (Fig^J, it is dif- 
ficult to determine the range of the orientation angle 9 
of these hydrogen-bonded species relative to the inter- 
face normal. The orientational distribution of these C2V 
species can be quite broad, different from that for the in- 
terfacial water molecules with the free 0-H bond. From 
our simulations, it appears to us that SFG measurement 
may not be very effective to determine the orientational 
angle of the C2V species at the air/water interface, even 
though it can do very well with the Coov 0-H bonds as 
shown above. However, our recent analysis of the SHG 
measurement of the neat air/water interface showed that 
SHG measurement might be able to help determine the 
orientational angle of the C2V species, but not the Cacv 



0-H bonds. Recent SHG results indicated that the aver- 
age orientation of the interfacial C2V water molecules is 
about 40° to 50° from the surface normaliS^ 

The molecular structure, orientation and dynamics 
at nonpolar material/water interfaces have been stud- 
ied by ab initio calculation, MD simulation, or them 
combined.^i^i^i^i^i^°i"i^'^i"i^°i^^ It appears that some dif- 
ferent conclusions were drawn on the molecular orien- 
tation and structure of the air/ water interface in dif- 
ferent studies&2iiS^ Nevertheless, many of these studies 
concluded that the dipole vector of the interfacial wa- 
ter molecules prefers lying parallel to the interface and 
have one of the 0-H bond protrude out of the liquid 
phase. The majority of the conclusions from theoreti- 
cal calculations agree satisfactorily with the experimen- 
tal analysis of ours and previous studies, but all the sim- 
ulation results were with significantly broader orienta- 
tional distributionsj'^1^1^^1^^1^^1^^1^'^1^^1^^1^'^1^^ There were 
reports concluded that some interfacial water molecules 
have their two 0-H bonds projecting into the vapor phase 
and with oxygen atoms in the liquid phase. How- 
ever, we have not found explicit spectroscopic evidence 
for such species at the air/water interface. These all in- 
dicate that detailed comparison of the theoretical cal- 
culations and the experimental analysis is certainly an 
important subject in the future studies. 



V. CONCLUSION 

Detailed understanding of the air/ water interface is 
important, and can be used for the general understand- 
ing of the liquid water structure. In this work, we pre- 
sented detailed analysis of the SFG vibrational spectra 
of the air/water interface taken in different polarizations 
and experimental configurations. Polarization and ex- 
perimental configuration analysis have provided detailed 
information on the orientation, structure and dynamics 
of the water molecules at the air/water interface. The 
success of these analysis indicated the effectiveness and 
ability of SFG-VS as a uniquely interface specific spec- 
troscopic probe of liquid interfaces and other molecular 
interfaces. It also indicates that for the neat air/water 
interface, as has been studied in the literature for some 
other simple air/liquid interfaces, the contribution from 
the interface region dominates the SFG spectra4£i2S42L2& 

Here are major conclusions we have reached for the 
air/water interface. Firstly, we concluded that the mo- 
tion of the interfacial water molecules can only be in 
a limited angular range, instead rapidly varying over a 
broad angular range in the vibrational relaxation time 
suggested previously. Secondly, because different vibra- 
tional modes of different molecular species at the in- 
terface has different symmetry properties, polarization 
and symmetry analysis of the SFG-VS spectral features 
can help assignment of the SFG-VS spectra peaks to 
different interfacial species. These analysis concluded 
that the narrow 3693cto~^ and broad 3550cm~^ peaks 
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belong to Coov symmetry, while the broad 3250cto~^ 
and 3450cm~-^ peaks belong to the symmetric stretching 
modes with C2v symmetry. Thus, the 3693cto~^ peak is 
assigned to the free OH, the 3550cm~^ peak is assigned 
to the single hydrogen bonded OH stretching mode, and 
the 3250cTO^^ and 3450cm~^ peaks are assigned to inter- 
facial water molecules as two hydrogen donors for hydro- 
gen bonding (with C2v symmetry), respectively. Thirdly, 
analysis of the SFG-VS spectra concluded that the singly 
hydrogen bonded water molecules at the air/water in- 
terface have their dipole vector direct almost parallel to 
the interface, and is with a very narrow orientational 
distribution. The doubly hydrogen bond donor water 
molecules have their dipole vector point away from the 
liquid phase. Finally, we did not find any observable ev- 
idence for interfacial water molecules with doubly free 
0-H bonds at the air/water interface. 

Many of the conclusions in this work agree well 
with previous reports, with much more detailed under- 
standings. The conclusion of the narrow range mo- 
tion of the free 0-H bond is different from the litera- 
ture. The explicit assignment of the broad 3550cm^^ 
peak and determination of the symmetry property of the 
hydrogen-bonded 0-H stretching modes in the 3100cto~^ 
to 3500cm^^ region are based on firm evidences. These 
conclusions as a whole provided a detailed and general 
picture of the spectroscopy, structure and dynamics of 
the air/water interface, which can be used for under- 
standing chemical and biological problems related to the 
ubiquitous water molecule in general. The concepts and 
approaches used in the analysis in this report can be ap- 
plied to studying on more complex molecular interfaces. 

Recently, extensive efforts with SFG-VS, as well as 
SHG, experimental studies and theoretical simulations 
have been devoted to the renewed interests on ion adsorp- 
tion and the Jones-Ray effect at the air/aqueous solution 

interfaceS.^4,64,74,99,lQ0,101,102,lQ3,104 g^ggggt ^^^^ 

tailed polarization and experimental configuration anal- 
ysis of the SFG vibrational spectra be applied to these 
interfaces. 
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Appendix: Calculation of d and c Parameters for 
Civ Molecule 

Here we present the expressions to calculate the pa- 
rameter c and d for water molecules with C21, symmetry 
using the bond polarizability derivative model first used 



by Hirose et aliS^iiS^ The detailed re-derivation of the 
complete expressions and the effectiveness of the model 
can be found in a recent review^ 

The relationship between the Raman depolarization 
ratio p and the bond polarizability r for a molecule group 
with G2V symmetry wasi^ 

^+^° (1-^)^(1 + 3 00=^ r) 

in which t is the H-O-H bond angle between the two 
OH bonds of a water molecule. With the Raman de- 
polarization ratio measured as about 0.03^2^ the bond 
polarizability r for OH bond in water molecule can be 
deduced to be 0.32, as used by Du et o/3 

The 7 hyperpolarizability tensor elements of water 
molecule with C2V symmetry are as the followings. 



Paac 


= ^"^^o^Ki + r) 


(1- 


r) cos t] cos( — 


Pbbc 


= ^^rcos( — ) 

UJal 2 






Pccc 


= ^'^^0^[(l+r) + 

UJal 


(1- 


r) cos t] cos( — 


Paca 


UJbl 


-r) 


sinr] sin( — ) 


Pbcb 









(13) 



Where G^ (1 + cosr)/Afo + 1/Mh and Gt = (1 - 
cosr)/Afo + ^/Mh are the inverse effective mass for the 
symmetric (ai) and asymmetric (61) normal modes, with 
Mq and Mff as the atomic mass of O and H atoms, re- 
spectively. UJal and ujhi are the vibrational frequencies 
of the respective modes. — 2^a^(^M^i as defined by 
Wei et alM' The water molecule are fixed in the molecu- 
lar coordination A'(a, 6, c) with the O atom at the coor- 
dination center, the molecule plane in ac plane, and the 
bisector from the oxygen to the two hydrogen atoms side 
is the c axis. 

For the achiral rotationally isotropic (Coou) liquid in- 
terface, the symmetric stretching (ss, ai) vibrational 
modes have;^ 

Axxz ■^yyz 

= ^^^(cos^ V')/3aac + (siu^ i)) (3bbc + /3ccc] (cos 6*) 

+ ^ A,[(sin2 ■4j)Paac + (cos^ i,)(3bbc - Pccc] (cos^ e) 

y.{2).ss ^ ^(2),ss ^ {2).ss ^ {2).ss 
Axzx ''^-zxx '^yzy '^zyy 

= -iA,[(cos0)-(cos3 0)] 

[(Sin^ ■Ip)f3aac + (COS^ tp)(3bbc - Pccc] 

xfh"' = [(sin2 ^)/3,,, + (cos2 i:)l3bbc] (cos 9) 

- Ns [(Sin^ ij)l3aac + (COS^ ^:)[ibbc - Pccc] {cos^ 9) 

(14) 
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And the asymmetric stretching (as, 61) vibrational modes 
have, 

{2),as ^ (2), as ^ (2),as ^ (2), as 
A.XZX Azxx -^yzy -^zyy 

- i A^s/3aca [(cos^ V) " (sin^ V)] (cos 9) 

+ iVs/3aca(sin2V')(cos3 6') 
^(2^as ^ 2iV,/3,,,(sin2V')[(cos0) - (cos^^)] (15) 

The 62 asymmetric mode are SFG inactive since the 
hyperpolarizabihty tensors fi^cb and (3cbb are zero. 

The Euler angel -0 can be integrated if the H-X-H plane 
of the XH2 group can rotate freely around its symmetry 
axis c. For water molecules with both OH bond hydrogen 
bonded to neighboring molecules in liquid phase, the Eu- 
ler angel V' should not be a fixed value. Assuming a ran- 
dom -0 distribution we have the following non-vanishing 
tensor elements for the symmetric-stretching mode.j^SiS 



but is very small. So the ss vibrational mode spectra in 
the sps and pss polarizations should vanish as the CH2 
group mentioned above. However, they have to be very 

TABLE III: The general orientational parameter c and the 
strength factor d for ss mode and as mode of water molecule 
with C2v symmetry in different polarization combinations. 
The d values bear the unit Pccc- 



ss mode 


d-ssp 


c-ssp 


d-sps 


c-sps 


d-ppp 


c-ppp 


Confif 




0.400 


0.038 


0.012 


1 


-0.146 


0.174 


Confif 


5-2 


0.374 


0.038 


0.013 


1 


-0.079 


0.338 


Confif 


5-3 


0.362 


0.038 


0.013 


1 


-0.046 


0.589 


Confif 


r.4 


0.257 


0.038 


0.012 


1 


0.066 


-0.378 


as mode 


d-ssp 


c-ssp 


d*c-sps 


c-sps 


d-ppp 


c-ppp 


Confif 


5-1 


-0.154 


1 


-0.122 


00 


0.262 


0.98 


Confif 


5-2 


-0.144 


1 


-0.129 


CX3 


0.272 


0.99 


Confif 


;.3 


-0.139 


1 


-0.128 


00 


0.277 


0.99 


Confie 


'A 


-0.099 


1 


-0.117 


00 


0.250 


1.01 



Xi't' = Xi'jr = -^N^iPaac + Pbbc + Wccc){cOs6) 

+ \N,{(5aac + Pbbc - 2[3ccc){cOS^0) 

v(2),^s = ^{2},'>s ^ (2),ss ^ (2),ss 
Axzx Azxx '^yzy -^^yy 

= -^N.iPaac + l%bc - 2l3ccc){{cOSe) - {cOS^O)) 
= lNs{Paac+Pbbc){cOsd) 

- ^N,iPaac+f3bbc-2p,,,){c0.3''e) (16) 

And the non-vanishing tensor elements for water 
asymmetric-stretching modes are, 

Xi'ir^ = X^^Jr = -lN,f3aca{{cOSe) - {cOS^ 9)) 

y.{2),as _ (2),as _ (2),as _ (2)ms _ Ij^ g / 3/,\ 
Axzx ~ /<.zxx ~ A-yzy Azya ~ 2 ' 

X^^Jr = N,l3aca{{cos9) - {cos'0)) (17) 

For CH2 group, there is a general relationship (3aac + 
Pbbc - 2/3ccc = 0, because r = This relation- 

ship makes x^x}x''=x^zxx''= xfzy''=X^zyy'' — 0' which 
means that the ss vibrational mode should vanish in the 
sps and pss polarizations according to Eq. 1161 For water 
molecule, r = 104.5°. Then hyperpolarizabihty tensors 
of the water molecule are as the foUowings: (3aac — 1.296; 

Pbbc = 0.557; Pccc = 1; Paca = Pcaa = 0.741; Pbcb = 

Pcbb = 0. Here all value are normalized to [3ccc — 1- 
Then, (iaac + Pbbc - ^Pccc = -0.147. This value is not 0, 



small comparing with in ssp spectra. This is fully con- 
sistent with the small intensities in the sps SFG spectra 
for the C2v water modes in Fig^ 

With above deduction, and following the proce- 
dure in previous report the general orientational 
parameter c and strength factor d for the symmet- 
ric stretching (ss) mode and asymmetric stretching 
(as) mode of water molecule in different polariza- 
tions and experimental configurations can be calculated 
(see Table 1111(1 . The parameters used in the calcula- 
tion are ni{Lu)—ni{uJi)—ni{LU2) = l; n2(ti;)=n2(wi) = 1.34; 
ri2(ijJ2)=l-18; n'{uj)=n'{LJi)— 1.15; n' {0)2)— 1-09, respec- 
tively. These parameters are the same as the dielectric 
constants used for calculation of the air/ water interface 
by Wei et al^ As we have discussed in our reportsj^S*^ 
polarization analysis with the co-propagating experimen- 
tal geometry is insensitive to the value of the dielectric 
constants of the IR frequency.^^'^'^ Therefore, we used the 
same refractive constants for the IR frequencies across 
the whole 3100cr7i^^ to 3800cm^^ region, and this does 
not appear to affect our analysis. These c and d val- 
ues are used to calculate the polarization and orientation 
dependence of the SFG intensity, as well as the inter- 
ference (phase) of different spectral features in different 
experimental configurations. These calculations can sat- 
isfactorily explain the detailed changes of the observed 
spectral features, as discussed in the main text. 

It is to be noticed that in the above discussion we only 
used single water molecule parameters. When there is 
association and clustering of water molecules, as long as 
the C2V symmetry preserves, and the H-O-H bond angle 
does not change significantly, above expressions dictated 
by symmetry properties should still be valid. 
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